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Determination  of  Kinetic  Parameters  of  Molecular  Beam 

Epitaxy 

Zhong  Jingchang 

(Changchun  College  of  Optics  &  Fine  Mechanics) 

Abstract 

A  kinetic  growth  model  for  molecular  beam  epitaxy  (MBE)  was 
discussed.  Furthermore,  high  energy  electron  diffraction  (HEED) 
was  used  as  a  surface  characterization  method  to  provide  evidence 
for  this  model. 

GaAs  was  used  as  an  example  to  study  the  growth  rate  of 
molecular  beam  epitaxy.  The  relation  between  the  growth  rate  and 
the  flux  of  Ga  was  verified  by  quadrupole  mass  spectroscopy. 

The  suitable  Ga  to  Al  partial  pressure  ratio  for  growing 
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materials  such  as  Ga(i_^l/^As  with  specific  values  of  x  was 
investigated  theoretically  as  well  as  experimentally. 

Furthermore,  the  relation  between  dopant  concentration  and  the 
corresponding  effusion  cell  temperature  was  given  for  Si,  Sn  and 
Be. 

X 

r 

I.  Introduction 

Molecular  bean  epitaxy  (MBE)  is  an  ultra  high  vacuum 
technique  to  produce  a  semiconducting  compound  through  the 
interaction  of  a  hot  molecular  beam  on  the  surface  of  a  crystal 
at  a  specific  temperature.  Unlike  liquid  phase  epitaxy,  MBE  is 
an  unequilibrated  growth  process.  This  special  mechanism  allows 
low  growth  rate  and  temperature.  Thus,  it  is  capable  of 
producing  large  area,  uniform  semiconductor  thin  films  (on  the 
atomic  level).  Furthermore,  the  thickness,  chemical  composition 
and  dopant  level  can  be  rigorously  controlled.  Important  MBE 
parameters  are  far  more  superior  to  those  of  liquid  and  gas  phase 
epitaxy  techniques  .  Because  of  this  reason,  it  is  capable  of 
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producing  structure  unattainable  by  the  latter  techniques. 
Therefore,  MBE  has  made  great  progress  in  the  fabrication  of 
microwave  components,  optoelectronic  devices,  integrated  optical 
devices,  and  the  theoretically  significant  quantum  trap  structure 
and  super  lattice  structure. 

The  MBE  growth  mechanism  is  being  studied  in  depth.  In  this 
work,  reaction  kinetics,  surface  morphology  and  MBE  growth  rate 
were  discussed  experimentally  and  theoretically.  Furthermore, 
the  relations  of  several  important  kinetic  parameters  were 
determined . 

II.  MBE  Growth  Model 


Advances  in  surface  analysis  created  favorable  conditions 

for  establishing  a  growth  model  for  MBE.  From  HEED  pictures  one 

can  observe  the  smoothening  process  on  the  substract  surface 

during  growth.  This  phenomenon  indicates  that  the  growth 

mechanism  is  realized  by  a  two-dimensional  stepwise  transport 

process.  As  atoms  adsorbed  on  the  surface  increase  to  the  brink 

of  forming  a  step,  any  point  on  the  surface  may  become  the  origin 

of  the  step.  The  repulsive  force  between  the  steps  creates 

acceleration  for  the  atoms  to  move  away  from  these  sources. 

Hence,  the  step  sources  migrate  on  the  surface  to  reduce  the 

[21 

density  of  steps  and  to  smoothen  the  surface 

Adsorbed  atoms  may  become  free  atoms  again.  This 
desorption  process  is  related  to  the  lifetime  of  adsorption  and 
the  concentration  of  adsorbed  atoms  according  to  the  following 
formula^ ^  : 
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where  p  is  the  desorption  rate,  n  is  the  concentration  of 
adsorbed  atoms,  and  t  is  the  surface  lifetime.  In  reality,  when 
a  beam  of  atoms  at  an  intensity  F  is  suddenly  injected  on  a 
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surface  in  MBE  the  rate  of  concentration  variation  of  adsorbed 
atom  is 


In  this  case,  the  initial  condition  is  n(0)  =  0.  If  we  assume 
that  r  £  f(n),  then  the  differential  equation  (2)  can  be 
integrated.  From  equation  (1),  we  get 


r-'t'-Fil-e'"*)' 

(3 

From  the  relation  between  t  and  temperature,  the  excitation 
energy  for  desorption  can  be  determined  by 


T-r0exp>  Ed  RT)  3 

It  was  experimentally  demonstrated  in  the  MBE  of  GaAs  that  the 
surface  lifetime  of  As  is  very  short  on  a  hot  clean  GaAs  surface. 
However,  it  is  relatively  long  on  a  Ga  covered  GaAs  surface.  The 
desorption  of  As,  in  this  case,  only  increases  with  the 
consumption  of  Ga.  This  indicates  that  it  is  necessary  to  allow 
Ga  adsorption  on  the  GaAs  surface  in  order  to  adsorb  As  at  high 
temperatures . 

In  the  actual  experiment,  elemental  arsenic  was  used  as  the 
source.  It  exists  in  the  AS4  form.  Molecular  beams  of  Ga  and  As 
react  on  hot  (100)  GaAs  surface.  AS4  molecules  were  first 
adsorbed  in  a  mobile  but  weakly  bound  state.  The  number  of 
adsorbed  Ga  atoms  controls  the  precipitation  and  reaction  of  As^. 
One  of  the  situations  is  that  it  is  desorbed  to  become  a  free 
As^.  The  other  is  to  be  chemisorbed  by  meeting  a  Ga  lattice.  Two 
neighboring  As^  molecules  in  the  Ga  lattice  may  react  in  pairs, 
i.e.,  the  surface  chemical  reaction,  to  form  a  new  desorbed  As^ 
molecule  and  four  dissociated  As  atoms  to  initiate  the 
chemisorption  reaction: 


’As,  tr'  — ►  As; -As*  D'  r, 
As:-MV\As>  4As‘;Ad), 


(B) 
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where  represents  the  surface.  If  an  adhesion  coefficient  - 
the  ratio  of  the  number  of  atoms  adhered  to  the  surface  to  the 
total  number  of  incident  atoms  -  is  used  to  describe  this 
process,  then  an  approximate  rule  may  be  obtained: 

whe"  r4S.»  2FGa'  1-e-'  SAs.<<0-5' 
then  GaAs  grows  on  the  surface  where  As  is  stable. 

When  FAs«  2FGa,  l.e.,  S^O.5, 
then  GaAs  grows  on  the  surface  where  Ga  is  stable. 

From  the  above  analysis  one  knows  that  the  rate  of  growth  of 

the  MBE  GaAs  thin  film  is  totally  determined  by  the  flux  F  of  the 

molecular  Ga  beam  incident  on  the  substrate  surface  (in 
2 

atoms/cm. s),  i.e.,  R  =  a.F  where  a  is  the  adhesion  coefficient. 

At  the  typical  MBE  growth  temperature  range  (450-620°C)  and  when 
As  is  stable,  its  numerical  value  is  close  to  1.  This  means  that 
almost  all  the  incident  Ga  atoms  are  combined  into  the  expitaxial 
layer. 

There  is  no  generalized  mechanism  for  the  incorporation  of 
dopants  in  MBE.  Therefore,  each  dopant  must  be  handled 
empirically.  Generally,  doping  is  realized  by  adding  the 
appropriate  element  in  an  effusion  cell  to  quantitatively  control 
the  electroactive  impurities  into  the  growing  thin  film.  The 
optoelectronic  property  of  the  thin  film  must  first  depend  on  the 
impurity  flux  across  the  growth  boundary,  the  probability  of 
atomic  adsorption  and  the  surface  lifetime.  Secondarily,  it  is 
also  related  to  the  actual  combining  behavior  and  the 
electroactivity  of  the  dopant. 

Sn  is  a  very  suitable  donor  dopant  for  GaAs  and  Ga^_xAlxAs. 

The  dopant  level  is  proportional  to  the  Sn  flux.  The  maximum 
19  -3 

level  is  10  cm  with  very  little  compensation.  However,  Sn  has 
the  tendency  to  precipitate  preferentially  on  the  surface.  The 
Sn  concentration  on  the  surface  is  several  orders  of  magnitude 
higher  than  that  in  the  bulk.  The  distribution  of  dopant  cannot 
vary  abruptly.  It  was  experimentally  discovered  that  the  donor 
level  increased  with  increasing  As-  flux  and  decreased  with 


rising  temperature  .  These  results  might  be  interpreted  by  the 
fact  that  the  combining  rate  for  Sn  is  related  to  the  Ga  vacancy 
concentration,  which  varies  with  the  As^  flux  and  the  growth  / 259 
temperature.  Sn  incorporation  is  also  limited  by  the  surface 
recombination  rate.  Before  a  steady  state  donor  concentration  is 
reached  in  the  GaAs  thin  film,  we  must  build  up  a  steady  state 
number  of  surface  Sn  atoms.  This  number  may  be  0.1  monolayer. 

Si  is  an  amphoteric  dopant.  In  GaAs  grown  by  MBE,  Si  is  a 

1 1 

donor  in  most  cases.  It  is  not  difficult  to  reach  a  flux  of  10 

2  16  18—3 

atom/cm  .s  to  obtain  a  doping  level  of  1  x  10  ^  5  x  10  cm"  . 

Be  is  an  ideal  acceptor.  In  the  incorporation,  there  are  no 
complications  such  as  preferential  precipitation  and  abnormal 
diffusion.  In  MBE,  the  free  acceptor  level  is  proportional  to 
the  temperature  of  the  Be  effusion  cell.  Therefore,  it  is  a  very 
promising  shallow  acceptor  in  the  GaAs  and  Ga^  xAlxAs  systems. 

III.  Experimental  Results  and  Discussion 

A  Perkin-Elmer  model  *  400  MBE  apparatus  was  used.  The 

10 

background  pressure  of  the  growth  chamber  was  less  than  10” 
torr.  An  ionization  tube  was  installed  near  the  substrate  to 
accurately  indicate  the  pressure  and  the  flux. 

1 .  The  Growth  Rate 

According  to  the  model,  when  the  temperature  remains 
unchanged  and  the  adhesion  coefficient  for  Ga  is  1  in  growing 
GaAs,  the  rate  of  growth  R  and  the  Ga  flux  F  have  the  following 
relationship : 

F-R>D-Nt/M-K'R,  R-T  t„ 
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Table  1  Growth  Rates  Calculated  Gallium  Fluxes,  and  Experi¬ 
mentally  Obtained  Currents  of  Gallium  Ionization 


1 .  substrate  number 

2.  growth  temperature  (*C) 

3.  partial  pressure  of  gallium  molecular  beam 

4.  partial  pressure  of  arsenic  molecular  beam 

5.  growth  thickness 

6.  growth  time 

7.  growth  rate 

8.  dopant  concentration 

9.  room  temperature  mobility 

10.  calculated  gallium  flux 

11.  gallium  current  measured  by  quadrupole  mass  spectrometer 


t  »  '  •  i  «  [  •  rm 
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where  D  is  the  density  of  GaAs,  Nq  is  the  Avogadro's  number,  M  is 
the  molecular  weight  of  GaAs,  T  is  the  thickness  of  GaAs  grown, 
and  tQ  is  the  growth  time.  Hence,  the  growth  rate  is  linear  with 
respect  to  the  gallium  flux.  Table  1  shows  the  growth  rates 
measured  at  487  +  3°C,  the  calculated  gallium  flux  values  and  the  /260 
ionization  currents  of  the  gallium  molecular  beams  as  measured  by 
the  Faraday  cylinder  in  the  quadrupole  mass  spectrometer.  The 
latter  two  should  only  differ  by  a  factor  determined  by  the 
instrument.  From  Figure  1  one  can  see  that  the  Ga  ionization 
current  is  also  linear  with  respect  to  the  growth  rate.  The 
difference  between  this  linear  realtionship  from  that  with  the 
calculated  Ga  flux  is  in  slope. 


< 
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Figure  1.  Dependence  of  the  Growth  Rate  on  the  Calculated 
Gallium  Flux  and  on  the  Experimentally  Obtained 
Gallium  Ionization  Current 


1. 

2. 

3. 

4. 


calculated  theoretical  Ga  flux 
measured  Ga  ionization  current 
Ga  ionization  current  (JO  '  A) 
Ga  flux  (10^  atom/s. cm^) 


5.  growth  ratio  (um/h) 
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Figure  2.  HEED  Photographs  Before  (left)  and  (right)  After 
Growth 

2.  Surface  Morphology 

Figure  2  shows  two  HEED  patterns  of  the  GaAs  (100)  surface 
in  the  MBE  growth  of  GaAs.  The  substrate  surface  was  finely 
polished.  However,  on  the  atomic  scale,  the  surface  is  still 
rough.  There  are  few  diffraction  fringes  and  they  are  coarse. 
After  2  minutes,  the  surface  became  smooth  on  the  atomic  scale. 
Very  fine  and  clear  diffraction  patterns  began  to  emerge.  This 
is  a  supporting  evidence  for  the  step  growth  and  surface  atom 
migration  model . 

3.  Ratio  of  Ga  to  Al  Partial  Pressure  in  Growing  Ga^  xAlx 

As . 

When  growing  Ga^_xAlxAs  with  a  specific  x  value,  the  ratio 
of  Ga  to  Al  partial  pressure  can  be  considered  as  follows: 

P  (;.»■•  _  L-.V  N  Ga'  r  M  <  M-  T  'r, 

P  .VI)  ”  A'  '  .v  Al  *.  )/  Al  '  I'  V  .  ’  (4) 

where  P(Ga)  and  P(A1)  are  the  partial  pressures  of  Ga  and  Al , 
respectively.  S(Ga)  and  S(A1)  are  the  sensitivities  of  the 


onization  gauge  to  Ga  and  Al,  respectively.  M(Ga)  and  M(Al)  are 
he  atomic  weights  of  Ga  and  Al ,  respectively.  T(Ga)  and  T(A1) 

.re  the  absolute  temperatures  of  the  Ga  and  Al  effusion  cells, 
espectively.  Table  2  shows  the  x  values  measured  by  electron 
:nergy  dispersion  (EDS)  spectroscopy  at  different  P(Ga)/P(Al) 
iartial  pressure  ratios.  If  the  x  values  in  Table  2  are  plugged 
nto  equation  (4),  by  taking  S(Ga)  =  1.7^®^,  S(A1)  =  0.9^^  and 
’(Ga)/T(A1)  =  1  to  represent  the  temperature  of  effusion  cells 
or  Ga  and  Al  in  the  MBE  process,  a  P(Ga)/P(Al)  value  can  be 
ibtained  corresponding  to  each  measured  x  value.  Table  3  shows 
.he  calculated  P(Ga)/P(Al)  values  from  equation  (4)  and  the 
elative  deviations  between  the  calculated  value  and  the  actual 
ralue  used.  One  can  see  that  the  values  given  by  equation  (4) 
ire  essentially  in  good  agreement  with  those  actually  used.  The 
lean  relative  deviation  is  3.2%. 

’able  2  X  Values  Measured  by  the  Electron  Energy  Dispersion 
Spectroscopy  in  the  Case  of  Different  Ratios  of 
P(Ga) / P ( A 1 ) 


s’  strate  temperature  (°C) 
G  partial  pressure  (Torr) 
A.  partial  pressure  (Torr) 
x  \alue  measured  by  EDS 


Table  3  Relative  Errors  Between  the  Calculated  P(Ga)/P(Al) 
and  Corresponding  Experimental  Values 
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1 .  x  value  measured  by  EDS 

2.  P(Ga)/P(Al)  used  in  experiment 

3.  P(Ga)/P(Al)  calculated  by  equation  (4) 

4.  relative  deviation 


Figure  3  Dependence  of  Si,  Sn,  and  Be  Doping  Concentrations 
on  the  Temperature  of  these  Dopants 

_3 

1.  dopant  level  (cm  )  _1 

2.  reciprocal  effusion  cell  temperature  1000/T(K“  ) 
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Figure  4.  Dependence  of  the  Experimentally  Obtained  Ga  and  Al 
Partial  Pressure  on  the  Temperature  of  Both  Elements 

1.  partial  pressures  of  Ga  and  Al  (Torr)  1 

2.  reciprocal  effusion  cell  temperature  (1000  x  °C“  ) 

4.  Relation  of  Si,  Sn  and  Be  Doping  Levels  and  Dopant  / 262 

Temperature 

Figure  3  shows  the  relation  between  the  concentrations  of 

Si,  Sn  and  Be  with  respect  to  the  effusion  cell  temperature.  One 

18—3  19 

can  see  that  the  level  could  reach  n  =  2  x  10  cm  for  Si,  n  =10 
cm"^  for  Sn  and  P  =  5  x  10^cm”^  for  Be. 

5.  Relation  Between  Partial  Pressures  of  Ga  and  Al  and 
Effusion  Cell  Temperature 

Figure  4  shows  the  relation  between  partial  pressures  of  Al 
and  Ga  with  respect  to  the  reciprocal  of  the  effusion  cell 
temperature.  This  curve  is  handy  in  establishing  the  required 
partial  pressure  ratio  by  controlling  the  temperature. 
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IV.  Conclusions 


In  molecular  beam  epitaxy,  the  crystal  growth  process  is 
kinetically  limited,  instead  of  quasi-equilibrium  limited  in 
liquid  phase  epitaxy.  At  a  fixed  temperature,  the  growth  rate 
exhibits  a  linear  relationship  with  respect  to  the  gallium  flux. 
This  experimental  result  is  in  good  agreement  with  the  theory  and 
serves  as  a  good  piece  of  evidence  for  the  growth  mechanism. 
Surface  characterization  by  HEED  and  quadrupole  mass  spectrometry 
also  have  significant  contribution  to  the  establishment  of  the 
model.  Although  only  GaAs  was  investigated,  it  might  be 
applicable  to  other  III-V  semiconductor  materials. 

Experimentally  measured  P(Ga)/P(Al)  ratios  agree  with  the 
projected  values  obtained  based  on  equation  (4).  From  a 
practical  angle,  this  is  a  good  reference  for  growing  Ga^_xAlxAs 
with  a  specific  x  value. 

This  work  was  conducted  under  the  direction  of  Dr.  M. 
Gershenzon.  Professor  of  Electrical  Engineering  and  Materials 
Science,  at  Southern  California  University  in  the  U.S.  Mr.  C. 
Shannon  also  assisted  in  this  work.  The  author  wishes  to  express 
his  gratitude. 
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